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We discuss the role of two different types of diagram in the proton-proton double parton 
scattering (DPS) cross section - single and double perturbative splitting graphs. Using 
T-H explicit calculations of simple graphs from these classes we show that the treatment of these 

graphs by the 'double PDF' framework for describing the DPS cross section, introduced 
a number of years ago by Snigirev and collaborators, is unsatisfactory. We suggest that 
^ ^ a contribution from single perturbative splitting graphs should be included in the DPS 

I cross section, albeit with a different geometrical prefactor to the contribution from 'zero 

perturbative splitting' graphs. 

1 'Double Perturbative Splitting' Diagrams in Double Par- 
^ ton Scattering 

We define double parton scattering (DPS) as the process in which two pairs of partons partici- 
pate in hard interactions in a single proton-proton (p-p) collision. DPS processes can constitute 
important backgrounds to Higgs and other interesting signals (see e.g. [T]), and can themselves 

fSj be considered as interesting signal processes, since they reveal information about parton pair 

correlations in the proton. 

Making the assumption that the hard processes A and B may be factorised, the cross section 

" . l for p-p DPS may be written as follows: 

X ^^(^,5) °^ X! / II c^2;aO?^b(Tjfc^A(s = a::ia;3s)o-j7^B(s = X2a;4s) (1) 

I a=l 

X ^ij{xi,X2,h; Q\, Q%)Tki{x3, X4, b; Q\, Qg) 

The cross section formula is somewhat similar to that used for single parton scattering 
(SPS), except that two parton- level cross sections a appear, and the PDF factors are two- 
parton generalised PDFs F (2pGPDs) rather than single PDFs. Note that in this formula the 
two 2pGPDs are integrated over a common parton pair transverse separation b. 

In many extant studies of DPS, it is assumed that the 2pGPD can be approximately fac- 
torised into a product of a longitudinal piece and a (typically flavour and scale independent) 
transverse piece: 

F,, (xi, X2, b; Qi, Ql) ~ D;^ixi,X2; Q%Ql)F{h) (2) 
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Figure 1: (a) A diagram that apparently contributes to the leading order DPS cross section 
according to the framework of [3]. The partons emerging from the grey proton blobs in the 
figure are nonperturbatively generated partons - i.e. ones existing at a low scale ~ ^qcd- (b) 
The 'crossed box' graph. In this part of the figure, A and B are arbitrary single particle final 
states with Q"^ > 



Then, if one introduces the quantity doff via (icff = !/[/ -F(b)^d^b], one finds that one may 
write cr^ entirely in terms of the longitudinal piece and a^s'- 

1 [ ^ 

'^iA,B)0^—^ I WdXaD''^{xi,X2;Q\,Q\)Dl\x3,Xi;Q\,Q\)cr^k-^Acrji^B (3) 

i.jjk.l a—1 

In [J a quantity D'^J {xi,X2]Q'^) is introduced, and an evolution equation for this quantity 
is given. We shall refer to the quantity and its evolution equation as the double PDF (dPDF) 
and the dDGLAP equation respectively. It is asserted in [3j that the dPDF is equal to the 
factorised longitudinal part of the 2pGPD in the case in which the two hard scales Q\ and 
are equal to a common value Q^. 

The dDGLAP equation contains two types of terms on the right hand side - 'independent 
branching' terms corresponding to emission of partons from a pre-existing pair, and 'single 
parton feed' terms corresponding to the perturbative generation of a pair from the splitting of 
a single parton. The single feed terms involve the leading twist single parton distributions as one 
might expect. Given this structure of the dDGLAP equation, with single feed terms included 
on the right hand side, a prediction of the 'dPDF framework' suggested in [3] for calculating 
the p-p DPS cross section is that a part of the 'double perturbative splitting' or 'Ivl' graph 
drawn in figure [l| a) should be included in the LO p-p DPS cross section. The part that should 
be included is proportional to [log((5^/A^)]"/(Teff at the cross section level, where A is some 
IR cutoff of order Aqcd, and n is equal to the total number of QCD branchings in figure [TJa) 
(including the two that only produce internal particles). This piece should be associated with 
the region of transverse momentum integration for the graph in which the transverse momenta 
of the branchings on either side of the 'hard processes' in the graph are strongly ordered up to 
scale ^ Q^. 

The question that then arises is whether such a structure in fact exists in the cross section 
expression for the loop of figure[l}^a). Starting from the conventional 'Feynman rules' expression 
for the loop, it is not immediately obvious what the answer to this question is. Here we 
will focus on answering this question for the specific very simple 'crossed box' loop shown 
in figure WPo) , which is predicted by the dPDF framework to contain a piece proportional to 
[log(Q^/A^)]^/crof^ The issues raised in the treatment of this example carry over to the more 
general loop of figure [I| a). 

^Note that the 'crossed box' topology of figure ifb) is the only box topology that contains a DPS singularity. 



2 



MPI@LHC 2011 



We expect the [log(Q^/A^)]^/(Tcff piece in figure[T]jb) to be predominantly contained in the 
portion of the cross section integration in which the external transverse momenta, as well as 
the transverse momenta and virtualities of the internal particles, are all small. This is actually 
the region around a certain pinch singularity in the loop integral known as the double parton 
scattering singularity |3] . In |S] , we obtained an analytic expression for the part of an arbitrary 
loop containing a DPS singularity associated with the loop particles emerging from the initial 
state particles being nearly on-shell and collinear, in the limit in which the external transverse 
momenta are small. Applied to the loop of figure [ijb) this reads (schematically, suppressing 
helicity and colour labels): 

Lopsfie^b) °^ / d^k<i>g^gg{x, k - Q2)$3^gg(l - X, -k) (4) 

X Mgg^A{s = X{1 - x)s)Mqq^B{s = x{l - x)s) + [q ^ q) 

In this formula, x = p2 ■ Qi/pi ■ P2, s — (pi +^2)^, and k (Q2) is the component of k 
{Q2) transverse to the axis defined by the directions of the incoming particles. $g_j.gg(a;, k) 
is the 0{as) hght cone wavefunction (LCWF) to produce a qq pair from a g [5], with the 
quark having lightcone momentum fraction x and transverse momentum k with respect to the 
parent gluon. It can be factored into a k and x dependent part, where the k dependent part 
is proportional to e • k/k^, e being the transverse part of the gluon polarisation vector. It is 
generally true that any QCD 1^-2 splitting with physically polarised external particles has a 
corresponding LCWF that is proportional to 1/k. Provided one uses a physical gauge for the 
gluon, all LCWFs corresponding to a 1 — )■ 2 QCD splitting are then proportional to 1/k. 

Inserting Q into the standard 2^-2 cross section expression, and performing a number 
of changes of variable, we arrive at the following expression for the DPS singular part of the 
gg — AB cross section: ^ 

CDPS.figgb) OC j Y\_ dXidXiaqq^A{s = XlXls)aqq^B{s = X2X2S) (5) 
-^V g^qq{xi,X2,r)T g^qq{xi,X2j -v) 



1=1 

r 



{2^f 

rg^qq{x^,X2,r) ^ -^6{1- X^ - X2)r^{x„X2) / d'k ]_^J ' j ' . (6) 

ZTT J [k+ 2 rj [k— 2 rj 

T''^ (xi, X2) contains a function of xi and X2 that may be regarded as a '1 — >■ 2' splitting 
function, multiplied by a constant matrix in transverse spac^ r is equal to the transverse 
momentum imbalance of one of the quarks/antiquarks in the loop between amplitude and 
conjugate, and is the Fourier conjugate variable of the parton pair separation b in the qq 
pair emerging from either gluon. Tg^qq{xi^ X2,v) can therefore be thought of as the 0(015) 
transverse momentum-space 2pGPD to find a qq pair inside a gluon. Note that the expression 
here effectively coincides with that of [7] , in which a cross section expression for the box of [T|b) 
is obtained starting from a pure DPS view of the box. 



^Note that the cross section is really a sum of terms with different T'^ {xi, X2) factors in the g qq 2pGPDs. 
This is associated with the fact that, from the point of view of the quarks, there is an unpolarised diagonal 
contribution to the process plus polarised and interference contributions in colour, spin, and flavour space. See 
e.g. [7l[8] for a discussion of correlation and interference effects in DPS processes. 
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Let us consider the part of the integral ([s]) that is associated with the magnitude of the 
imbalance r being smaller than some small cut-off A that is of the order of Aqcd- The con- 
tribution to the cross section from this portion contains a log^(<5^/A^) factor multiplied by 
(which can be thought of as an effective 'l/cTcff' factor for this contribution). The majority of 
this contribution comes from the region in which the transverse momenta and virtualities of 
the quarks and antiquarks in the gg — >■ AB loop are much smaller in magnitude than \/Q^ (i.e. 
the region in which the assumptions used to derive Q apply), which is a necessary feature of 
a contribution to be able to regard it as a DPS-type contribution. By making a specific choice 
of A (let us call this As), one could obtain an expression which is exactly in accord with the 
expectations of [3j - that is, a product of two large DGLAP logarithms multiplied by the same 
l/ccff factor that appears in diagrams in which the parton pair from neither proton has arisen 
as a result of one parton perturbatively splitting into two ('2v2' or 'zero perturbative splitting' 
diagrams). The 1/acs factor for the 2v2 diagrams presumably has a natural value of the order 
of that is set by the nonperturbative dynamics {Rp = proton radius). 

The fact that we have to make a somewhat arbitrary choice for A in order to arrive at the 
result anticipated by the dPDF framework is concerning. There is nothing in the calculation 
of figure [ijb) to indicate that we should take the region of it with |r| < As as the 'DPS part' 
- the scale As does not naturally appear at any stage of the calculation. There is no more 
justification for taking the part of the box with |r| < As to be the DPS part than there is for, 
say, taking the piece with |r| < 2A5, or that with |r| < A5/2, to be the DPS part. We have 
had to artificially introduce the cut-off A5 in figure [T](b) in order to obtain a power-suppressed 
DPS part because there is no scale in the graph apart from Q^, so in order to obtain a term 
proportional to l/Q^, a second scale has to be introduced 'by handj^ 

There therefore appear to be some unsatisfactory features of the dPDF framework with 
regards to its treatment of the box in figure [IJb). In a physical gauge, precisely the same issues 
will be encountered for the case of the arbitrary 'Ivl' graph in figure [T|a). One obtains a result 
that is consistent with the dPDF framework if one demarcates the portion of the cross section 
integral in which the transverse loop momentum imbalance between amplitude and conjugate 
is less than As as DPS, but there is no natural reason to do this. There is no distinct piece of 
figure [ija) that contains a natural scale of order Aqc d and is associated with the transverse 
momenta inside the loop being strongly ordered on either side of the diagram. In fact, most 
of the contribution to the total cross section expression for the graph comes from the region 
of integration in which the transverse momenta of particles inside the loop are of 0(^7 Q^l^ 
This fact suggests that at the level of total cross sections, we should perhaps remove 'Ivl' 
graphs entirely from the DPS contribution, and regard them as pure SPS (this approach is 
advocated in p2j, for both the total and the differential cross sections). Treating the graphs 
in this way would have the advantage that we would not perform any double counting between 
DPS and SPS - the graph of figure [l}^a) is in principle also included in the SPS pp — > AB cross 
section (albeit as a very high order correction that will not be included in practical low order 
calculations, if the number of QCD emissions from inside the loop of the graph is large). 

Very similar conclusions may be reached if one uses a covariant gauge such as the Feynman 
gauge for the gluon fields in figure [IJa), although these conclusions are perhaps not obtained 

•^This is related to the fact that in massless perturbation theory, there are no power corrections. 

*One should bear in mind however that the same is not true for the cross section expression differential in 
the transverse momenta of A and B for Pj- y^^Pj^ b Here, the major contribution is associated with 

transverse momenta in the loop <C \/ if there is one emission or more from inside the loop, or with a range 
of transverse momenta between ~ |pt,a|, |pt,s| and ~ \fQ^ if there are no such emissions [5]. 
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Figure 2: (a) A generic graph from the '2vl' class, (b) A simple '2vl' graph in which a gluon 
splits into a qq pair, and these partons then interact in two separate Drell-Yan interactions with 
a 'nonperturbatively generated' qq pair from a proton. 

so readily. In a covariant gauge, gluons with unphysical 'scalar' polarisation can exist in loop 
diagrams. Such scalar-polarised gluons can give rise to power-law DPS divergences rather than 
logarithmic ones, and additional 'super-leading' contributions to the AB production process 
(in terms of powers of Q) - the two phenomena are related. On the other hand one generally 
expects the 'super- leading' contribution to cancel in a suitable sum over graphs (as in [9]), 
which effectively leaves one with the same logarithmic DPS divergences that are encountered 
in a physical gauge. 

One can gain some insight into the source of the problems in the dPDF framework by 
looking at the b-space 2pGPD corresponding to T g^qq{xi,X2,h). This comes out as being 
proportional to - this behaviour (which was first spotted in [7]) can be traced to the fact 
that the g qq LCWF in b space (like any any LCWF corresponding to a QCD 1 — )■ 2 splitting 
with physically polarised external particles) is proportional to 1/b, and r(b) ~ 4>(b)^. Note 
that this behaviour is very different from the behaviour of all 2pGPDs that is anticipated by the 
dPDF framework (i.e. smooth function of size Rp). There is no natural feature in the product 
of two 'perturbative splitting' 2pGPDs that is of size Rp and can be naturally identified as DPS. 
A key error then in the formulation of the dPDF framework seems to be the assumption that 
all 2pGPDs can be approximately factorised into dPDFs and smooth transverse functions of 
size Rp. A sound theoretical framework for describing p-p DPS needs to carefully take account 
of the different b dependence of pairs of partons emerging from perturbative splittings, whilst 
simultaneously avoiding double counting between SPS and DPS. 

2 'Single Perturbative Splitting' Diagrams 

Aside from the 'Ivl' graphs that were the focus of the previous section, and the '2v2' graphs 
that were also briefly mentioned, there is a further class of graph that can potentially contribute 
to the p-p DPS cross section. These are graphs in which one proton provides one parton to 
the double scattering, and the other two, at the nonperturbative level - a representative graph 
is sketched in figure [2ja). For obvious reasons, we will refer to the graphs as '2vl' or 'single 
perturbative splitting' graphs. 

It seems clear that we should include contributions from the 2v2 graphs as part of the DPS 
cross section. An important question is whether we should also include contributions from the 
2vl graphs, and if so, what form these contributions should take (in particular, how does the 
effective (Toff factor for the 2vl graphs differ from that for the 2v2 graphs?). 

To answer this question, let us take a similar approach as we did for the Ivl graphs in the 
previous section. That is, we take the graph drawn in figure [2](b) that has the simplest possible 
2vl structure, and see whether there is a 'natural' part of the cross section expression for it 
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that is proportional to l/Rp, and also contains a large logarithm associated with the 1 — 2 
splitting. If there is such a structure in the 2vl graph, then this part of this graph should be 
included in the LO DPS cross section, and we also expect there to be a \og{Q^ / A^)" / piece 
in the more general 2vl diagram of figure ^a) that should also be included in the LO DPS 
cross section. 

In the calculation of the cross section for figure [2f^b) , it is necessary to include a hadronic 
amplitude or wavefunction factor ipp to find two nonperturbatively generated ('independent') 
partons in the proton, at the amplitude level in the calculation. It would be inappropriate to 
try and calculate a 2vl cross section in a naive 'fully parton-level' way omitting the proton on 
the 'nonperturbative pair' side because then one has three particles in the initial state (whereas 
the standard framework for calculating a cross section requires two particles in the initial state) . 
Furthermore, by deleting the proton on the 'nonperturbative pair' side one is then neglecting 
the important fact that the nonperturbatively generated partons are tied together in the same 
proton (as was pointed out in |10|). The use of proton wavefunctions or hadronic amplitudes 
in the calculation of DPS-type graphs was discussed long ago in [TT], and has been discussed 
more recently in ^ [TU] . 

After a lengthy calculation, one finds that the cross section for figure contains the 
following expression: 

a'li,2(s) =&qq^.y*is = Xiyis)&qq^^^{s = X2y2s) (7) 

7(2^ ^''"^^'^^''^^''"^] ^^3^99(2^2)0(1-2/1 - 1/2) y^^ —2 



In this expression we have omitted helicity and colour labels and sums for simplicity. The 
quantity Tp-qq {xi, X2','r) is the 2pGPD of the nonperturbatively generated parton pair, whilst 
Pg^qq IS thc LO 1 — >■ 2 Splitting function for the process g — >■ qq. m is a symmetry factor that 
is equal to 1 if the two hard processes are the same, and is equal to 2 otherwise (for the double 
Drell-Yan process under consideration, it equals 1). 

The integral over Qi in Q gives rise to a large transverse momentum logarithm log{Q^ / A'^) , 
whilst the integral over r supplies a prefactor of order ^ 1 /R^ (we assume ^p-qq (r) only 
has support for |r| values of order Aqcd - see later). Thus, there is a part of the cross section 
expression for figure[2](b) that is proportional to \og{Q^ / A^) / and should be included in the 
LO DPS cross section. 

When we generalise the result ([t]) to the leading logarithmic part of the arbitrary 2vl 
diagram in figure [2ja), and then sum up all of the diagrams to obtain the contribution of 2vl 
graphs to the LO DPS cross section, then we obtain the result belov 



o-^_B) (s) =2 X y / dxidx2dyidy2&tk^A{s = xiyis)aji^Bis = a;2y2s) (8) 



^2j^^2^p!i"dep(yi; 2/2, T; Q^) 



The quantity D'^J {xi,X2\ Q^) is the 'accumulated sPDF feed' contribution to the dPDF. This 
evolves from a zero initial value at a nonperturbative scale Qo ^ Aqcd according to the full 

^Note that for simplicity we take the two hard scales to be equal here, = = Q^, and only write down 
the unpolarised diagonal contribution in colour, flavour and spin space. The contributions associated with spin 
polarisation (longitudinal or transverse) are expected to have a similar structure. On the other hand, it is known 
that the colour and quark number interference contributions will contain Sudakov logarithms - see e.g. (71 15). 
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dDGLAP equation. Tp\^j^^p{yi,y2,r; Q^) is the 'independent branching' 2pGPD. This sums up 
the effect of independent strongly ordered parton emissions from a nonperturbatively generated 
parton pair. It evolves according to the dDGLAP equation with the sPDF feed term removed. 
There is an additional prefactor of 2 in (|8]) because there are two sets of lv2 graphs that give 
equivalent contributions ~ in one set the nonperturbatively generated parton pair emerges from 
the 'left' proton, whilst in the other it emerges from the 'right' proton. 

A critical requirement for the derivations of Q and ^ to be valid is that parton pairs 
connected only via nonperturbative interactions should have an r-space distribution that is 
cut off at values of order Aqcd (or equivalently a b-space distribution that is smooth on 
scales of size <^ Rp)- This appears to be a somewhat reasonable requirement - at a low scale 
Qo ^ A.QCD there is only the scale Aqcd available to set the size of the r profile for Tp\ndep^ 
and the evolution of Tp^indep essentially preserves the transverse profile to higher scales. What 
is more, such behaviour for Tp i^dep would appear to be required in order to get the necessary 
prefactor of order 1/Rp in the 2v2 contribution to DPS, which is calculated according to the 
following expression (for the diagonal unpolarised contribution): 

'^iA^B)(^) J dxidx2dyidy2^ik^Ais = xiyis)&ji^B{s = X2y2s) (9) 



X J (l#r;^,»<iep(^i-^2,r;Q^)r^:„,,,(yi,2/2,-r;g^) 

Note that the quantity {2n)-^ J rf^rr^J^„^^p(yi, 7/2, r; Q^) in ^ is equal to fp'i„dep(yi, y2, b = 
0; Q^). This appears to indicate that the 2vl contribution to DPS probes independent branching 
2pGPDs at zero parton separation. In fact, the result ([s]) corresponds to a broad logarithmic 
integral over values of b^ that are <C but » The quantity f ^'-^^^^(yi, ?/2, b = 0; Q^) 

appears as a result of our smoothness assumption on rp'.i„dep(2/i' 2^2, b; Q^). 

If one assumes that the independent branching 2pGPD can be approximately factorised 
according to the prescription in ([2]), then the contribution to the DPS cross section from 2vl 
graphs is similar to that predicted by the dPDF framework, albeit with a different associated 
'cTeff'- Indeed we find that (cre//,2t)2)"^ — J d'^h[F {h)]"^ , whilst (cre//,it,2)^^ ~ F{h = 0). If one 
then makes the further assumption that F{h) is approximately Gaussian, one finds that each 
2vl contribution to DPS receives a factor of 2 enhancement over the 2v2 contribution from 
the (cToff)^^ geometrical prefactor (as is also noted in [10]). One should perhaps not put too 
much trust into this exact figure, however - it clearly relies on a number of assumptions whose 
validity is somewhat uncertain. 

3 The Total Cross Section for DPS 



If one were to take the suggestions outlined earlier in this report at face value, then one would 
obtain the following expression for (the unpolarised diagonal contribution to) the total LO DPS 
cross section: 



(^(A,B){s) = cr, ; (s) + a,,g. (s) (10) 



with cr^^^^(s) and cr^^^^(s) being given by the expressions ([8| and ^ respectively. This 
expression shares some common terms with the DPS cross section formulae proposed in [lOj 



and |12j . Looking closely at (10) however, one can identify a number of concerning issues with 
regard to this equation, which indicate that modifications to it may be required in order to 
correctly describe the DPS cross section. 
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The first issue is that we were originally expecting to obtain an expression for the DPS 
cross section looking something like ([T]), with the 2pGPDs in these formulae each having an 
interpretation in terms of hadronic operator matrix elements. Our proposed expression ( 10 ) 
deviates somewhat in structure from these expectations (in particular, one would struggle to 
come up with a matrix element representation for D^J (xi,X2]Q^))- This feature is related to 
the fact that we have entirely removed the 'Ivl' contribution from the DPS cross section. 

The second issue is that there is a rather sharp distinction in ( 10 1 between perturbatively 
and nonperturbatively generated parton pairs, with the 2pGPD for the latter Pp'-^^^p having a 
natural width in r space of order A. Does there exist some scale at which we can (approximately) 
regard all parton pairs in the proton as being 'nonperturbatively generated' in this sense (as is 
assumed in (10))? If so, what is the appropriate value for the scale (presumably it should be 
rather close to Aqcd)? 

A final issue is that in the above, we have largely ignored the interesting and potentially 
important interference and correlated parton contributions to DPS catalogued in [71 [5] ■ 

In this report, we have shown that the treatment of Ivl and 2vl contributions to DPS 
by the dPDF framework of Snigirev et al. appears to be unsatisfactory, and presented the 
results of a calculation that indicates that we should include a contribution to the DPS cross 
section from 2vl graphs, if we include a contribution from 2v2 graphs. There appear to be 
some unsatisfactory features in our alternative suggestion for the DPS cross section (10), which 
perhaps indicates that completely removing any contribution from Ivl graphs from the DPS 
cross section is not quite the correct prescription. 
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